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Abstract 

We examine the contributions to rare processes that arise in models where the 
Standard Model fermions are localized at distinct points in compact extra dimensions. 
Tree-level flavor changing neutral current interactions for the Kaluza-Klein (KK) gauge 
field excitations are induced in such models, and hence strong constraints are thought 
to exist on the size of the additional dimensions. We find a general parameterization 
of the model which does not depend on any specific fermion geography and show that 
typical values of the parameters can reproduce the fermion hierarchy pattern. Using 
this parameterization, we reexamine the contributions to neutral meson mixing, rare 
meson decays, and single top-quark production in e + e~ collisions. We find that it is 
possible to evade the stringent bounds for natural regions of the parameters, while 
retaining finite separations between the fermion fields and without introducing a new 
hierarchy. The resulting limits on the size of the compact dimension can be as low as 
TeV" 1 . 
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1 Introduction 



In recent years there has been much interest in the possibility that there may exist compact 
extra dimensions with sizes far above the Planck length. In particular, the possibility of 
TeV _1 -sized extra dimension arises in braneworld theories PU El El lU El E] ■ By themselves, 
they do not allow for a reformulation of the hierarchy problem, but they may be incorporated 
into a larger structure in which this problem is solved, such as the case of large extra 
dimensions El E]- I n the scenario with TeV -1 extra dimensions, the Standard Model 
(SM) fields are phenomenologically allowed to propagate in the bulk. These models are hence 
subject to stronger experimental constraints and have distinct experimental signatures from 
the case where gravity alone is in the bulk. 

There are many possibilities for how to place the Standard Model fields in the TeV -1 
bulk. In the universal extra dimensions scenario all fields see the extra dimensions, giving rise 
to a conserved parity that relaxes direct production and precision electroweak constraints, 
and may provide a dark matter candidate [TTH ITT) IT^l ITUj . The effects of universal extra 
dimensions in rare processes have been considered in f3J E3 EH QZj- It is also possible 
to localize the fermions without localizing the bosons, which allows for the gauge fields 
to propagate freely throughout the bulk. More recently it was noticed by Arkani-Hamed 
and Schmaltz (AS) that one could localize different fermion species at different points in 
the TeV -1 extra dimensions ^H]- These "split fermion" models naturally suppress many 
dangerous operators, particularly those inducing proton decay. They also can naturally 
generate large Yukawa hierarchies; and it has been shown by multiple authors that there 
exist models which can generate the correct spectrum of fermion masses, as well as the 
correct magnitudes for CKM matrix elements [Till EHl EU E21 EH El] ■ The most stringent 
generic limits in this case arise from precision electroweak measurements, which place the 
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compactification radius at R < 2 — 4 TeV -1 [23 I2E1 123 121] • The specific fermion locations 
can be probed in high energy collisions, and at very large energies, cross sections will rapidly 
vanish since split fermions will completely miss each other in the extra dimensions. |29t 13*0] 

This makes the split fermion scenario an attractive possibility for the origin of the 
Yukawa hierarchy However, split fermions (like most models of the fermion spectrum) are 
also capable of generating large flavor changing neutral currents (FCNC). The magnitude 
of these currents in the neutral meson sector has been estimated by several groups, and 
apparently generate strong constraints |2H IBTj. I32j . In this paper we reexamine these compu- 
tations to derive more model independent constraints on split fermion models arising from 
FCNC and show that it is possible to evade the stringent bounds for natural regions of the 
parameters. 

This paper is organized as follows. In Section 2 we set up the split fermion scenario 
in as much generality as possible and give statistical arguments to demonstrate that they 
can account for the observed fermion spectrum. We then describe how FCNC are generated 
in this scenario. In section 3 we calculate the effects on neutral meson oscillation. Section 
4 presents the effects on rare B decays, and single top production in e + e~ collisions, and 
Section 5 concludes. 

2 The Model of Split Fermions 

Here, we construct a very general model that characterizes the effects of separating the 
Standard Model fermions in an extra dimension. We start by examining the original model 
considered by AS. 

In the AS model there is one extra dimension, which is taken to be flat. It is possible 
that this extra dimension is actually a "brane" with a finite width embedded in some other 
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extra-dimensional scenario. For this reason TeV -1 dimensions are often called "fat branes", 
but they need not be tied to other models. Note that if the brane is not a string theory 
object, but arises from some field theory mechanism, then it necessarily has finite extent 
in the extra dimensions. This makes the study of fat branes essential to building realistic 
field-theoretic models of extra dimensional scenarios. 

In this model the Standard Model fields are localized to the brane. Note that the 
word brane here refers to any mechanism for achieving this localization. It may or may not 

be the same as the branes encountered in string theory. Initially, all fields are allowed to 

propagate in the entire dimension. In addition to the fields present in the Standard Model, 

there is a real scalar field which couples to the fermions, but not to the gauge bosons or the 

Higgs. 

If the scalar has a ^-symmetric potential, then it can develop a stable solution which 
tunnels from one of the vacua to the other, called a kink solution. A mechanism for localizing 
fermions to a thin but finite width region inside a domain wall has been known for some time 
[33j. There it was noted that in 1+1 dimensions a massless fermion with a Yukawa coupling 
to a scalar field that has a kink-profile vacuum expectation will develop a zero mode, with a 
Gaussian profile centered at the location of the kink. This can be trivially extended to more 
dimensions by considering a domain wall instead of a soliton and making all zero modes 
constant in the transverse directions. Note that a five-dimensional fermion field contains 
two four-dimensional fermions, one of each chirality. If the extra dimension is infinite, then 

the zero-mode of only one chirality is normalizable. If the extra dimension is finite then 
something else is needed to produce chirality. A standard procedure is to compactify the 
dimension on an S l / Z2 orbifold, which projects out the unwanted chirality. A nice side-effect 

of this is to render the kink absolutely stabled 

Ut is interesting to think that if one invokes a mechanism to localize the gauge bosons, as in then 
one could have a fat brane residing in an infinite dimension. 
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In contrast to the fermion sector, the gauge bosons are free to propagate throughout 
the extra dimension. Since the dimension is compact, and flat, the mass spectrum of the 
Kaluza-Klein gauge states is linear with M 2 = n 2 /R 2 , and the orbifold boundary conditions 
project out the odd solutions, so the wavefunctions along the fifth dimension, y where < 
y < R, are 



where R is the size of the extra dimension. Putting all this together allows investigation of 
brane world models where there is a single extra-dimension of roughly inverse TeV size with 
fermions localized in the center and gauge bosons propagating though the entirety. 

A more interesting picture can be obtained by thinking about the fermion localization 
mechanism. There is a simple heuristic for why this should occur. The fermion is Yukawa 
coupled to a scalar field which develops a non-zero VEV. The ordinary fermion Higgs phe- 
nomena should then give the fermion a mass. However, the VEV is position dependent and 
in particular there is a place where it is zero (the center of the kink). So the fermion has 
a position dependent mass, which is somewhere zero. Thus, the fermion is easiest to excite 
near the zero mass, and so most of the probability for the lowest lying state (the zero mode) 
will live near the center of the kink. 

Given that heuristic, it should be reasonable that if the 5D fermion has a mass M, then 
the center of the Gaussian moves toy = M/2/z 2 , where fi is the slope of the kink profile, and v 
is the scale of the VEV. Indeed, it turns out that this is the case, as was first noted by Arkani- 
Hamed and Schmaltz [18 . This allows different fermion fields to be localized at different 
points in the extra dimension. To see why this is desirable, consider an operator, O, that 
involves fermions separated by a distance d. The effective 4D coupling in the dimensionally 
reduced theory is proportional to the integral over the extra dimensions of the wavefunctions 
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of all fields appearing in O. Since the fermion wavef unctions are Gaussian, this gives a 

suppression proportional to e _ct/i , where a depends on the operator being considered. This 
has been shown to be very effective at suppressing dangerous higher-dimensional operators, 
such as proton decay. Additionally, the fact that exponentially different couplings can result 
from linear separations provides a natural means of explaining the fermion mass hierarchy. 
Lighter fermions have greater separation between their left and right handed components. 
In this way Arkani-Hamed and Schmaltz proposed a theory to explain the Yukawa hierarchy 
without invoking new symmetries, and which is safe from proton decay. Several authors 
have proposed specific "geographies" that do indeed reproduce the correct fermion masses, 
as well as the CKM parameters [151 12H1 12T] . 

There are, however, other potentially dangerous effects of the fermion separation 

which are not suppressed by this mechanism. The gauge bosons will have a Kaluza-Klein 

(KK) tower of states. The zero modes, which are flat in the extra dimensions, correspond 

to the SM gauge fields, and have the correct couplings to the fermion zero modes. On 

the other hand, the excited states have cosine profiles, as given in Eq. JU).* The coupl me; 

strength of these modes to the fermions are scaled by an integral over the overlap of the 

fermion and gauge wavef unctions. However, since the height of the boson wavefunction will 

be different at the locations of the different fermions, there will be non-universal couplings 

of a single gauge KK-state to different fermion species. This leads to the possibility of flavor 

changing interactions, including tree-level neutral currents, for the KK-modes of the 7, Z, 

and gluon, as illustrated in Fig^ One then expects large effects to come from the tree- level 

contributions of the KK gluon states to FCNC processes, in particular to neutral meson 

oscillation. Calculation of these effects can put limits on the size of the extra dimension. 

Also, note that while this discussion was motivated by the kink model, these issues will be 
tin general they are exponentials, e my , but the orbifolding projects out the odd modes. 
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relevant to any model with split fermions. This is an example of the general principle that 
any attempt to explain the Yukawa hierarchy will necessarily treat flavors differently, and 
will tend to generate large flavor- changing effects. 




Figure 1: Feynman diagram for the tree-level KK gauge exchange mediating neutral meson oscil- 
lation. 

In practice, geography independent constraints have been difficult to obtain due to 
the large number of parameters in the model. These are, R, a, the width of the fermion 
wavefunctions (which is 1/// in the kink model), and (d n — d) positions, where d is the number 
of extra dimensions and n is the number of independent fermion fields. Previous discussions 
PH EH E2] have put constraints on R only by first obtaining a single set of positions that 
reproduce the Yukawa couplings of the Standard Model, and calculating the flavor-changing 
effects in that particular geography. However, one would like a more model-independent way 
of understanding the magnitude of flavor effects in this class of models. 

To accomplish this we consider the problem of FCNC in split fermion models in as 
much generality as possible. A specific, realistic model exists in string theory [32], as well as 
the field theory example just presented. In summary, we abstract from these the following 
points: 
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1. There exist one or more extra dimensions, compactified with a radius R.§ 

2. Each fermion field, tpi has a chiral zero mode that is localized near the center of the 
dimension at yj, with Gaussian profiles ip ~ e ~(y-yi) 2 / a2 ^ where the width a <C R. 
If there is more than one extra dimension they are taken to be isotropic in those 
dimensions ip ~ e -(y-yi) 2 /°\ 

3. The gauge bosons are free to propagate in the entire "fat brane" part of the extra 
dimension. (There may be a larger bulk accessible to gravity.) 

4. The boundary conditions for the bosons are taken to be such that the wavefunction 
for the n-th KK-mode is A ~ cos(nny/R). Note that these are generally the same 
conditions that allow chiral zero modes for the fermions. 

5. The field content (gauge group, number and charge of matter fields) is identical to the 
Standard Model, plus whatever fields are necessary to localize the fermions. 

These assumptions generate an effective four-dimensional Lagrangian that reduces to the 
Standard Model at low energies. The new features present are the propagating gauge KK- 
modes, their couplings, and the fact that the Higgs Yukawa couplings are determined by the 
fermion locations. 

We now construct the interaction Lagrangian for this scenario, focusing on the quark 
sector in this paper. An analogous treatment of the leptonic sector can be performed. 
Note that there are excited states of the fermion fields in addition to the KK boson states. 
However, since the fermions are localized with a width smaller than R, the scale of the fermion 

excitations will be significantly higher than that of the KK gauge states. In addition, the 

fermion KK modes do not participate in the processes considered here. We therefore only 

§In one dimension the compactification is S 1 /Z2. In more dimensions we take the compactification to be 
fiat, and orbifolded in such a way that it looks like a simple product of single dimensions. 
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consider the fermion zero modes, while we include the complete KK-tower for the bosons. 
With one extra dimension the coupling of the n-th KK boson to a flavor localized at the 
scaled position i = x/R is determined by the overlap of wavefunctions 

/ dy^(y)ij(y)A {n \y)~ [ dy cos (nny) e -(^W- 2 ~ cos ( ravr £) e ~n^ 2 /R\ ( 2 ) 
Jo Jo 

where y has now been normalized to R. For 5 extra TeV -1 dimensions this generalizes to 

l ' R \ (3) 



c^(£)= (j[coa(n k ir£ k )j e~ r ^ 2 



The gauge coupling of the gluons, for instance, can then be written in flavor space as 

£ int = V2g a Gf n)li (d L ^T A C^d L + d R ^T A C^d R ) + (d -> u) + h.c. (4) 

Here &l(r) is the vector of left (right) handed down-type quarks, d = (d s b), g s is the SU(3) 
coupling constant, and G^ 11 is the n-th KK gluon field. The diagonal matrices are the 
wavefunction overlaps given by Eq. ©. The factor of y/2 arises from the rescaling of the 
gauge kinetic terms to the canonically normalized value for all n. 

Now, the Higgs zero mode, which is the Standard Model Higgs, is flat in the extra 
dimension, H° ocl/R. Then the Yukawa couplings to the 4D Higgs field are given by 

RX 5 [ dyH q L q R ~\ 5 [ dy e -^-Vi) 2 l°\-(.v~V3? ~ \ e ~^- y ^ la \ (5) 
Jo Jo 

Here A5 is an overall 5D coupling constant that is fixed to be 0(1) by the top quark mass. 

We write the 4D Yukawa couplings to (for instance) the down-type quarks in the flavor basis 
as 

^Yukawa = dV^ M d V^ d (6) 

Where MdVjf is the matrix of Yukawa couplings with elements given by Eq. (|5|). and 
Md is the diagonal mass matrix. 



We can now write the relevant terms of the Lagrangian as 

£ = d L V^M d vl d) d R + U L V^M U V[ U) "R + ^W^u L rd L 

oo 

+ [^9sG^ )A (d L YT A CPd L + d R ^T A C^d R ) + (d - u)] + h.c. (7) 

n=l 

After the usual transformation to the mass basis, the CKM-matrix is clearly the product 
y^y( d \ Note, however, the presence of non- universal couplings prevents the products 
jj<i{ n ) = ytenQWyW f rom De i n g trivial, so there are flavor- changing interactions in the KK- 
gluon sector. These also occur in the excited photon and Z couplings. However, those are 
suppressed relative to the gluons by a factor of g/g s , so we expect that the KK-gluons will 
dominate any process to which they contribute. 

Before examining the numerical impact of the tree-level FCNC interaction in rare 
processes, it will first be useful to get a handle on how far the fermions need to be separated. 
It has been shown by Grossman and Perez JH] that there exists at least one set of positions 
that correctly reproduces the observed fermion spectrum and magnitude of the CKM ele- 
ments. They found that, subject to a certain set of naturalness assumptions, there was a 
single solution. A different solution was found in [21] by choosing different up and down-type 
Yukawa coupling constants in the 5D theory. Typical separations in these solutions are from 

1 — 20 units of the fermion width. In what follows, we parameterize the separation between 

2 fermions in units of the width, i.e. Ay = y% — yj = a^a, and treat as phenomenological 
parameters. In addition, we find it useful to define p — cr/R. 

As a counterpoint to the studies in and [21] we have performed a simple Monte- 
Carlo analysis in an attempt to see how large of a hierarchy is generated naturally for fermions 
randomly distributed on an interval. To do this we randomly draw fermion positions from 
a distribution flat on the interval [0,a max ], and use these to compute the Yukawa matrices 
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from Eq. (jSJ). We then compute the singular values of these matrices, which are the fermion 
masses. We can get a sense of the hierarchy by taking a particular Yukawa matrix (say the 
up-type) and finding the ratio of the largest to smallest singular value. In Fig |2] we show a 
histogram of the log of this ratio for ct max = 15. For comparison we have computed the same 
value for a "null hypothesis" where instead of the split fermion scenario, the entries of the 
Yukawa matrices themselves are drawn directly from a distribution flat on the interval [0, 1]. 
As expected, the case of split fermions clearly generates a much larger hierarchy. What is 
surprising is that one needs to set a max ~ 10 — 15 before a hierarchy of six orders of magnitude 
becomes common, while in it was claimed this hierarchy could result from a max ~ 5. 
The discrepancy is due to the fact that, while a separation of a = 5 will indeed generate 
a matrix element of order 10 -6 , the singular values (which are the actual masses) of a full 
Yukawa matrix with separations no larger than 5 will tend to be too large. We note that the 
full fermion spectrum can be generated by p as large as 1/15, i.e. without introducing a new 
large hierarchy between the compactification and fermion localization scales. Also note that 
a max represents the part of the extra dimension in which the fermions can be localized and 
need not be the same as 1/p, which is the size of the dimension through which the gauge 
bosons can propagate. 

3 Constraints from Neutral Meson Oscillation 

Significant effects from the flavor-changing gluonic couplings should show up in neutral 
meson oscillations. We start by examining the effects on Kaon mixing. The AS* = 2 effective 
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Log 10 (m 1 /m 3 ) 

Figure 2: Histograms of 10 5 random trials of the size of the fermion hierarchy, log(mi/?rt3), where 
mi is the largest and the smallest mass for a single Yukawa matrix. Main graph: In the split 
fermion model where the positions are randomly drawn on the interval [0, 15] in units of the fermion 
width. Inset: In a "null hypothesis" where the Yukawa matrix elements are randomly drawn from 
the interval [0, 1]. 



Lagrangian from the single KK-gluon exchange depicted in Fig. ^ is 

9 °° 1 

^eff — n9 8 / j M 7, / , U i(sd) U j{ds) a *l S * a 3<» S 3 



r s ^ mi 

n=l n i,j=L,R 



i,j=L,R 



E(cos(mrx di ) - cos(rmx Si )) (cos(rmx d .) - cos(mrx s ,)) 2 
m 6 i • (8) 

n=l ™ 
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Here, the X{ are the positions of the d quark fields, yi of the s quark, and we have used the 
unitarity of the Vs (here, we have approximated this with 2x2 unitarity, we discuss the 
third-generation effects below). We are especially interested in the form of the sum over the 
KK-modes. While KK-sums are usually divergent in more than one extra dimension and 
require a cutoff, ours contains a natural cutoff arising from the finite width of the fermion 
zero mode, and hence is convergent for a number of extra dimensions 5 > 1. This is simple 
to understand physically. The cutoff sets in when the wavelength of the KK-mode is of order 
the fermion width, which occurs at -R/n max = a. At higher momenta the wavefunction of 
the boson is oscillating many times within the fermion allowing it to resolve the fermion's 
wavefunction, and exponentially decouples. The fact that this cutoff arises naturally in 
the field theory model is an attractive feature of that particular mechanism for fermion 
localization. 

In one additional dimension the sum converges even without the exponential sup- 
pression. In this case it is insensitive to the value of p and can be computed analytically by 
ignoring the exponential factor. To do this we need to evaluate 

F(x,y)^£ (C ° S(ri7rx) 7 OS(ri7ry))2 , (9) 

71=1 



and 



G(x y) - ( cos ( n7ra: i) ~ cos(mryi))(cos(mnK 2 ) - cos(mn/ 2 )) ^ 



n 

n=l 



The computations for these sums are presented in the appendix. The final result is 

F{?,y) = —\x-y\, (11) 



and 
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G(x 1 ,y 1 ,x 2 ,y 2 ) = yOi ~ X A + \Vi -Jfel - ki -Jfel - \%2 -yi\)- (12) 
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This tells us that the flavor changing effects depend, as expected, on any nonzero separation 
between fermion fields. 

The hadronic matrix elements for the gluonic contributions to Kaon mixing are given 
by (computed in the vacuum insertion approximation) 



(K°\d L Ys L d Ll ,s L \K°) = X -f K m k (13) 
{K°\d Ll »s L d R lMK Q ) = f 2 K m k (± + \ 



as well as those with (L <-> R), which have the same evaluation. Written out in full, the 
contribution to Attik is then 

Am K =ne(K°\£ AS=2 \K ) 

= ^g 2 R 2 ( \Vi n Vt 12 \ 2 F{x dL ,x SL ){K Q \d L ^sdLlM^ 

+ |Vg u V£ 12 \ 2 F(x SR ,x SR )(K°\d R Ys R d Rl ,s R \K°) (14) 
+ {Vt uV^yf nVg 12 )G(x dL , x SL , x dR , x s R )(K°\ d L *f8i4w R \K°) 

+ {V£ lx V R * l2 Vt\ 1 V£ l2 )G{x dL ,x SL ,x dR ,x s ^ 

where x dh R are the positions of the d field, and the x SL R of the s. Note that all possible sep- 
arations (between quark fields) are present, but some enter with different signs. In principle 
then, the gluonic contribution could be made small for any values of R and p by placing 
the quarks at appropriate places. However, the terms involving only right or left handed 
fields occur with the same sign. So, to achieve significant reduction, cancellations must occur 
between those terms and the terms which involve both chiralities. This involves tuning the 
quark positions to the values of different hadronic matrix elements, which introduces a fine- 
tuning problem. Otherwise, it would imply that the UV physics that localizes the quarks has 
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information about the IR behavior of QCD! We can therefore expect that cancellations will 
be 0(1) at most. This is seen clearly in the Monte Carlo trials, where the random positions 
have no relation to the hadronic matrix elements and no significant cancellation occurs. 

In light of this we can explore the magnitude of the flavor effects just by looking 
at a single term in (jl4j) ; for convenience we choose the first. We can then describe the 
contributions with only three parameters: the radius R, the scale ratio p, and the separation 
between one pair of fermions, a. The sum over KK- modes is then calculable in terms of 
these parameters. Since R enters only in the mass in the KK propagator, we can write the 
contribution in the simple form 



where V ds ds stands for the appropriate product of 4 elements of the V^' R matrices, two at 
each vertex, and, as shown in F(x,y) only depends on the difference of it's arguments, 
so we can write it as a function of only a single variable, given by the product pa. The 
subscript on F reminds us that in two dimensions or more F also depends on p directly as 
the cutoff parameter, in which case it must be computed numerically. If we demand that 
this contribution to Amx be no larger than the measured value (a conservative assumption 
from the point of view of constraining the model) we get 



where (3k is a coefficient of dimension 1, which depends on the meson parameters. This 
expression immediately generalizes to other neutral meson systems by using the appropriate 
coefficient (3, and the appropriate matrix elements of V. Table Q shows the values of /3 for 
cases of interest, along with representative values of F p (pa). 

The resulting constraints are shown in Fig. |3] and EJ for 1 and 2 extra dimensions 
respectively, using the value of f3 and V£ s ds appropriate for the Kaon sector, and assuming 



A™x = -g 2 J K m K R 2 Vl As F p (pa) 



(15) 




(16) 
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a=y/Rp 

Figure 3: The behavior of the constraint on 1/R with a for various values of p. The area below 
the curves is excluded. Note that the size of the additional dimension in units of the fermion width 
is 1/p, so the curve for p = 1/10 ends at a maximal separation of a = 10. 

that the V are CKM-like in magnitude (we discuss that assumption in detail below). There 
are two features of note. First, with one extra dimension the constraint is a simple square- 
root function, as can be seen from Eq. (jllj) . This means that the flavor- changing effects 
can be made arbitrarily small by reducing p. That is, by increasing the hierarchy between 
fermion and boson scales. Second, in two dimensions the effect seems to be roughly constant 
in p, and flattens off at large a. We know that the sum over the KK states is diverging 
logarithmically before it reaches the cutoff, and so it should be getting larger as p decreases. 
However, shrinking p brings the fermions closer together making the flavor effects smaller. 
In two dimensions these two effects are seen to roughly cancel. In three or more dimensions, 
the divergence of the sum wins completely, and the bounds on are huge, effectively 
removing these cases from consideration as realistic models. 

In Figs. 03 and |H1 we display our results for all meson mass differences, taking p = 
1/100. For mixing in the Kaon sector and B® sector, the bound is set by demanding that the 
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Figure 4: Constraints on 1/R from Kaon mixing for two extra dimensions; the area below the 
curves is excluded. 

new physics produce an effect which is no larger than the observed value. For D° mixing, 
the effect is restricted to lie below the current experimental bound. There is no experimental 
upper bound on B° s mixing, so we assume two values, one the size expected in the Standard 
Model, the other about 4 times larger, corresponding to the curves labeled small and large, 
respectively. Note that the most stringent constraints come from mixings involving the first 
and second generation. 

This pattern suggests a loop-hole in the otherwise stringent constraints. Namely, the 

V matrices need not be CKM-like. Since the CKM matrix is the product V^V^ d \ the 
observed CKM hierarchical structure could result from a completely different structure at 
the level of the V 1 . If there is small first to second generation mixing the Kaon and D° 
constraints will be relaxed, and all constraints would then be of order a few TeV, even for 
large values of p. 

However, in the previous calculation we ignored the third generation when imposing 
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Figure 5: Constraints from all species of neutral meson mixings for one extra dimension, taking 
p = 1/100. See text for a description of the experimental values used. Note that the K° and D° 
results are separate lines that overlap due to a numerical coincidence. 
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Figure 6: Neutral meson constraints for two extra dimensions. The area below the curves is 
excluded. 
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p 


a = 1 


a = 5 


a = 10 


a = 20 


a = 50 


ID 












10- 1 


1.15 


2.33 


4.76 


NA 


NA 


io- 2 


0.036 


0.23 


0.48 


1.00 


2.45 


10~ 3 


0.0012 


0.023 


0.048 


0.097 


0.25 


2D 












io- 1 


2.05 


3.82 


3.59 


NA 


NA 


io- 2 


2.25 


5.27 


6.46 


7.47 


8.22 


10~ 3 


2.27 


5.41 


6.67 


8.08 


9.78 



Meson 


(3{ MeV) 


K° 


1125.86 


B° d 


478.01 


B° s 


67.4346 


D° 


1124.23 



Table 1: Left: Representative values of the sum F p (pa) for one and two extra di- 
mensions. Right: Multiplicative j3 factors for mass splittings of the neutral mesons; 
l/R>Py/V*F p {pa). 

the unitarity condition in Eq. (jHJ). Transitions between two 4D mass eigenstates will involve 

all three generations in the localization (flavor) basis. In this case the matrices Jj[ will 
contain the positions of all three generations of quarks, and the unitarity conditions on the 
V matrices will be changed. For instance, for the term in Eq. (|T4*|) with both left-handed 
chiralities (and dropping the L index), in place of \Vj* u V^* 12 \ 2 F(x d , x s ) we should have 

, X s , Xfj, X s ) 

+V 3 \V 32 V 1 * 3 V n G(x b ,x s ,x d ,x s ) + \V 31 \ 2 \V 32 \ 2 F(x b} x s ) (17) 

These additional terms (including the ones not displayed above corresponding to right and 
mixed chiralities) will insure that in any mass splitting observable many mixing angles and 
fermion separations will enter. It then becomes non-trivial to reduce Am^ by adjusting 
mixings alone. However, it is still possible to reduce the Kaon and D° constraints by noticing 
that the new terms contain fewer diagonal elements. Hence, if the weak and mass eigenbases 
are not too badly misaligned, i.e. the V 1 have large diagonal elements and smaller off- 
diagonal elements, then the strongest constraints may be relaxed somewhat. To get a better 



sense of what is typically possible, we again run Monte Carlo simulations. From these we 
learn that a typical suppression factor is 10 _1 for the factors multiplying (3 in Eq. ([16)1 . 10~ 2 
is not uncommon, and 10 -4 is obtainable, but rare, even for the fairly large value a max = 15. 
This is illustrated in Fig. 
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Figure 7: Histogram of values of the factor H = v/J V*F(x, y) summed over the appropriate 
fermion positions and V matrix elements. 



All of these considerations together show that once the parameter space is thoroughly 
explored, it is possible to evade the large constraints from meson mixing for natural regions 
of the parameters. 



4 Rare Decays 

We also consider processes involving only a single flavor- changing vertex. The best examples 
of this type which receive contributions from KK gluon exchange are rare B decays, such 
as B — > ipK$ and B — > <pKs- The most interesting aspects of these decays are, of course, 
their associated CP-violating asymmetries. However, since we have no control of the phases 
present in split fermion models, we can't address the new contributions to CP-violating 
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observables in a model independent fashion. Nonetheless, there are tree-level strong coupling 
contributions to these decays, so we can expect significant contributions to the branching 
fractions in split fermion models. 

The effective Lagrangian for a process with a single flavor change is given by 

9 °° 1 

£es =1 = ig 2 s jp E uSva%w%- (is) 

n =l n i,j=L,R 

Since one vertex is flavor diagonal, the sum depends on the absolute position of the fermions. 
For instance, the analog of Eq3. (f9llTT|) is 

cos(n7rz)(cos(n7rx) — cos(nTry)) 



00 



F\x,y,z) = J2 



n=l 
r 2 



n 2 



— ( \z + x\ + \z — x\ — \z + y\ — \z — y\ + ttx 2 — -ny 2 ) . (19) 



where z corresponds to the location of the quark at the flavor conserving vertex. This 
additional complication turns out to be minor, as the actual magnitude of the sum is similar 
to that in the previous case and does not vary much over the parameter space, as can be 
seen from Fig. |HJ 

We consider the decay amplitude 

A(B^<pK s ) = \g 2 s R 2 ^ F'{x bi , x Si , x Sj )(^|(^7^)( Wi)|£°> (20) 

* i,j=L,R 
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Figure 8: KK sum for processes involving a single flavor changing vertex as a function of one of the 
fermion positions, corresponding to the case where one of the fermions is localized at the orbifold 
fixed point. The position of the third (flavor conserved) fermion, z, is varied in steps of one unit, 
with z = corresponding to the top curve and z = 20 to the bottom. The sum is done for a single 
extra dimension of size 20 units. 

The relevant matrix elements are jHHl 

(<bK s \(s L rb L )(s Llll s L )\B°) 

(4>K s \(s L Yb L )(s R ^s R )\B°) 
((t>K s \{s L ^b L , 3 ){s R ^,s R ^\B Q ) 

Here, i,j are color indices, displayed explicitly in the non-singlet terms which now contribute. 
The common factor is 

H = 2(ei-p B )f+mlF + (ml), (22) 

where represents the polarization vector of the phi meson, and F + (q 2 ) is the form factor 
for this decay. There are an additional four matrix elements obtained by taking (L <-> R). 



-H 



H 



12 



(21) 
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We use the values / = 233 MeV [37|, F + (7nJ) = 0.38 [HE], and m = 1020 MeV [3H]. 

For the branching fraction we obtain (ignoring the 0(1) differences among matrix 
elements) 

/ OO Jjd] rjd \ 2 

and similarly 

^-^^S.GxlO-^^Vp^)^ J] i^*, . (24) 

\i,j=L,R / 

Demanding that this not be larger than the observed rate gives the approximate constraints 
1/R > 1.0 TeV from B — > <f)K$ and 1/R > 0.5 TeV from B — > ipK$. These are not 
competitive with those from B^ and B s meson oscillation, and provide a good consistency 
check. It is interesting to note that if a way can be found to reduce the constraints from the 
Kaon sector to the few TeV scale without disturbing the 6-quark couplings (say by arranging 
the mixing angles), then this contribution to B — > tpKs,4>Ks is of roughly the same order 
as that of the Standard Model. Any new phases in this scenario will thus contribute to the 
CP-violating observables with equal effects in each decay channel. 

We have also estimated the ree-level KK contribution to single top quark production 
at LEP, e + e~ — > iq + tq (with q = u,c) which proceeds via KK r y^ rl ' and exchange. Using 
the parameterization in jlU] we find an effective anomalous flavor-changing vector coupling 



22 



vz given in this case by 

OO TTU Tjl 

v ZcS =(V2 + 2 sin 9 W Q) M^B} y £ -J^LJM (25 ) 

i,j=L,R n=l 

^6.9xl0- 4 f - X y (^$4)1 y F'(x c ,x t .,x e .)] , (26) 

where we have used only the c — > t since it dominates the u — > £ by a factor of 10. The 
current constraint on this effective flavor changing coupling from LEP data is vz < 0.75 |41j . 
and so split fermions in the up-quark sector are not significantly constrained. 

Lastly, we examine the potential gauge KK contributions to the purely leptonic decay 
of neutral B mesons, B® — > £ + £~, with q = d,s. This process proceeds through tree- level 

7^ n ) and KK exchange, and could ,in principle, probes fermion splitting in the leptonic 
sector as well as in the quark sector. Here, for simplicity, we will assume that the leptonic 
fields are all localized at the same point in the extra dimension and will ignore any possible 
flavor changing leptonic interactions. The Ab = 1 Lagragian describing this decay is then 

c?£= 1 = y 2G l £ w E u$$)fatfhi glJuA , (27) 

a=y,Z n=l n i=L,R 

where G 1 = e and Gz = g/ cos8 w . The hadronic matrix element governing this decay is 
given by 

<o|<?! )7/ X (n W^ + ^7/X {n W^°> = ifMg b L , 7/ zU d L {n) - &, 7/z u$ n) ] , (28) 

where represents the momentum of the B° meson. Due to parity, only the axial-vector 
current contributes to this matrix element. In the case of photon KK exchange, such con- 
tributions are generated when the left- and right-handed fermions are localized at separate 
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points. Here, we will ignore this possibility and consider the case where only the ex- 
change mediates this decay. The branching fraction is then 



B{B S -> n+ff 



7 r B K 

7TCf„ 



1 - 



Ami 



m 



B 



1/2 



n- 



1.15 x 10" 



R 



T/4 

db,fifi 



lTeV" 1 ; (0.04) 2 



9L,zF'(x bL ,x qL ,xe L ) — gR,zF'(x bR ,x qR ,X£ R ) 



(29) 



where we have chosen g = s and used = 200 MeV. The function F' is as given in Eq. (19) 
and has a value in the range —2 to +2 as shown in Fig. 8. Taking values of this function 
which maximizes the sum over the KK states, i.e., F' = 2, yields a value of unity for the 
sum, resulting in a branching fraction of 1.15 x 10~ 6 for R = ITeV -1 . This is a significant 
enhancement over the Standard Model value[42 j of B(B S — > /U + /i~) — 4.0 x 10 -9 . The 
experimental bound on this decay, B{B S — > /i + /U~) < 2.6 x 10~ 6 , as determined by CDF[43j, 

sets the limit > 815 GeV when the sum over the KK states takes on its maximal value. 
The sensitivity of this decay mode in probing the size of the additional dimension is thus 
comparable to that of B s mixing and will improve with Run II data at the Tevatron[~ 



5 Conclusions 

Models where the Standard Model fermions are localized at specific points along a compact 
extra dimension offer an attractive means for constructing the fermion mass hierarchy and 
suppressing dangerous operators such as proton decay. In these scenarios, the fermions obtain 
narrow Gaussian wavefunctions in the additional dimension with a width much smaller than 
the compactification scale. The fermion Yukawa couplings are then generated by the overlap 
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of the localized wavefunctions for the left- and right-handed fermions. Lighter fermions are 
thus more widely separated than heavier ones. 

The gauge fields are free to propagate throughout the bulk in these scenarios and 
their KK excitations develop tree-level flavor changing interactions which are proportional 
to the overlap of their wavefunctions with those of the localized fermions. Gluons, as well 
as the electroweak gauge bosons, then mediate flavor changing neutral current processes at 
dangerous levels. Previously, it was thought that the only way to avoid stringent bounds 
on the size of the compact dimensions was to minimize the separation of the fermion fields, 
thus endangering the scenario's natural explanation of the fermion hierarchy. 

In this paper, we have reinvestigated these new FCNC interactions and have per- 
formed a general, systematic, model independent analysis. Our results hold for any such 
model of the fermion hierarchy with specific fermion geographies. We have employed a model 
parameterization which contains only three parameters: the size of the extra dimension R, 
the scaled width of the localized fermion p = a/R, and the fermion separation distance 
expressed in units of the width, Ax = acr. We performed a simple Monte Carlo analysis and 
determined that the fermion mass hierarchy can be reproduced in our parameterization for 
natural values of the parameters. 

We then evaluated the KK gluon tree-level flavor changing contributions to neutral 
meson oscillations. We found that the sum over the KK states is exponentially damped 
for higher KK gauge states as the KK states can then resolve the finite size of the fermion 
wavefunction. This allows us to perform the KK sum in a scenario with more than one extra 
dimension. We then evaluated the constraints from Kaon mixing in the case of one extra 
dimension and confirmed previous results that l/R ^ 100's TeV for larger values of p unless 
the separation was very small. However, the constraint shrinks to l/R > few TeV for smaller 
values of p, even for widely separated fermions, at the expense of introducing a new hierarchy 
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between the compactification and fermion localization scales. The constraints from Bd and 

B s mixing were found to be much less restrictive. We also performed the evaluation for the 

case of two or more additional dimensions and found that the FCNC constraints were much 
more difficult to evade. 

We next studied the dependence of our constraints on the fermion mass mixing matri- 
ces, and found that with a realignment of the matrix elements our bounds could be reduced 
further by factors of 10-100. 

In addition, we examined the rare meson decays Bd — > ipKs, 4>Ks, as well as single 
top-quark production in e + e~ collisions, and found weaker limits of the size of the extra 
dimension of order TeV -1 . We note that the KK gluon contributions to these rare decays 
are significant enough to generate interesting effects in the related CP violation observables. 

In summary, we have shown that once the parameter space is systematically explored, 
it is possible to evade the stringent bounds from FCNC in split fermion models for natural 
values of the parameters and without the introduction of any additional hierarchies. Lastly, 
we note that the introduction of brane localized kinetic terms are known to significantly 
reduce the couplings of gauge KK states |H3 EH] and may help to even further reduce the 
constraints from FCNC in these scenarios. 
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6 Appendix 



Here we present a cute way to perform the sum over KK modes analytically in one dimension, 
and see that the sum is exactly linearly proportional to the separation. [17] The functions 
that we need are 

(cos(n7rx) — cos(mn/)) 2 

n=l 



F{x,y) = J2 



rr 



(30) 



and 



G(xi,x 2 ,yi,y 2 ) = E 



n=l 



^cos(n7ra;i) — cos(n7n/i))(cos(wr22) — cos(n7n/ 2 )) 

n 2 



(31) 



We can do both of these by evaluating 



/(*) = E 



cos m 



n=l 



rr 



(32) 



So that 



F(x, y) = ~(/(27ra) + f(2ny) + y - 2f(nx + iry) - 2/(ttx - iry)) 



(33) 



G(x 1 ,x 2 ,yi,y 2 ) 



f('KX 1 + nx 2 ) + f(nxi - nx 2 ) + f(ny 1 + ny 2 ) + f(-Ky 1 - iry 2 ) 



- f(irx 1 + ny 2 ) - f(-Kx 1 - iry 2 ) - f{-Ky 1 + ixx 2 ) - f(ny 1 - 7rx 2 ) 



(34) 



Writing the cos (no;) as two exponentials and combining the sums we have 



2" 



(35) 
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The function / is then the solution of the differential equation 

1 00 1 

rw = -^^ = -- emx + o ( 36 ) 



2 2 

n=£0 n=— 00 



-i (2tt f; *(a;-27rA0-lJ • 



(37) 



This is solved by 



H(x) = -- (n\x\ - — - — ) , (3J 



2 V 2 3 



where H is defined on the interval [— ir, 71] and is 27r-periodic for other values (this takes care 
of the sum over delta functions). We then have 

f(x) = H(x) + ax + p. (39) 

Looking at the original function we see that we must have f_ f = which gives (3 = 0. 
Also, since / is an even function a = 0. 

We can then use Eq. ()39|) in and and use the physical condition that all 
arguments are positive to get the final result 

F{x,y) = —\x-y\, (40) 



and 

r 
4 



G{xi,x 2 ,yi,y2) = -— ( \%i - x 2 \ + \y% - y 2 \ - \%i - 2/2 1 - \yi — 2/2 1 } - (41) 
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